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We study the linearized Bogoliubov excitation spectrum of infinitely long anharmonically trapped 
Bose-Einstein condensates, with the aim of overcoming inhomogeneous broadening. We compare 
the Bogoliubov spectrum of a harmonic trap with that of a theoretical flat-bottom trap and find 
a dramatic reduction in the inhomogeneous broadening of the lineshape of Bogoliubov excitations. 
While the Bragg excitation spectrum for a condensate in a harmonic trap supports a number of 
radial modes, the flat trap is found to signiflcantly support just one mode. We also study the 
excitation spectrum of realistic anharmonic traps with potentials of finite power dependence on the 
radial coordinate. We observe a correlation between the number of radial modes and the number 
of bound states in the effective potential of the quasi-particles. Finally we compare a full numerical 
Gross-Pitaevskii simulation of a finite-length condensate to our model of infinite, linearized Gross- 
Pitaevskii excitations. We conclude that our model captures the essential physics. 



I. INTRODUCTION 



The spectrum of weak excitations over the ground state of a Bose-Einstein condensate (BEG) has been found to 
obey the Bogoliubov dispersion relation 0. There is however, a finite width to the dynamical structure factor S{k,u>) 
which characterizes the response of the BEG to an exciting field with momentum k and frequency to. Usually the 
exciting field is in a Bragg configuration 2], where k is the momentum difference of two laser beams, and uj is the 
detuning between the two. There is a small intrinsic width originating from the Beliav decay of quasi-particles which 
is usually dominated by broadening due to finite time of the Bragg pulse and inhomogeneity of the condensate. The 
latter heavily depends upon the trap geometry. In some cases the wavelength of an excitation is much shorter than the 
condensate dimensions. Under these circumstances it is justified to employ the local density approximation (LDA) 0. 
For cigar shaped condensates, the radial dimension is sometimes comparable to the excitation wavelength, leading to 
the breakdown of the LDA, and the appearance of radial modes 0,t-5]. Avoiding inhomogeneous broadening is possible 
in the spectral domain using echo spectroscopy where a degenerate Bragg pulse transfers Bogoliubov excitations with 
wavenumber +k to —k 0, and in the time domain where rapid oscillations cause a suppression of inhomogeneous 
mean field and of Doppler dephasing 0. Suppression of the inhomogeneity mechanisms, allows longer coherence times, 
and opens the possibility of studying the homogeneous broadening mechanisms, that reflect the intrinsic decoherence 
processes of the bulk excitations, e.g. elastic collisions with the BEG @. 

In this paper we propose the use of a flat-bottom trap in order to reduce the inhomogeneous broadening of Bogoli- 
ubov excitations in a BEG. We restrict ourselves to the cigar shaped geometry, in which the inhomogeneity is mainly 
manifested by the presence of many radial modes. By linearizing the Gross-Pitaevskii equation (GPE), assuming a 
uniform and infinite potential in the axial direction and using a quasiparticle projection method .5], we obtain the 
Bogoliubov excitation spectrum for a flat-bottom trap. Then, we compare the above results to the excitation spectrum 
of a harmonically trapped BEG. A clear reduction in the inhomogeneous broadening is evident by comparing between 
the two spectra. This result, although intuitive, is not trivial for two reasons. The first is that even in a flat trap the 
density is not constant, but varies over a distance characterized by the healing length ^. 

The second reason is that even a hypothetical homogeneous finite condensate supports several radial modes. Indeed, 
we find that it is the gradual decay of the density to zero at the trap boundary that reduces the coupling to the high 
order radial modes. The envelope of both harmonic and flat trap spectra are found to match the LDA We then 
examine the excitation spectrum of a condensate trapped in the power-law potential V — Hpff , where p is the radial 
coordinate, and p is an even integer. We observe the substantial excitation of a second radial mode in the spectrum, 
and demonstrate its suppression using higher order potentials. The transition between potentials supporting one or 
two radial modes is explained by an effective potential picture. Finally, we solve the full GPE and find the Bogoliubov 
spectrum for a realistic, experimentally feasible Q three-dimensional axially elongated trap, confined both in the 
radial and axial direction using a high-order potential. The spectrum obtained by the linearized approach where the 
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potential is axially uniform, and by the full GPE for an elongated condensate, are found to be qualitatively similar 
due to the tight radial confinement. 

The outline of the paper is as follows: in section we study the ground state and the quasiparticles of the flat- 
bottom trap, then in scction llTTl a comparison between the flat-bottom trap and the harmonic trap excitation spectra 
is given. In section Hvl wc study high order anharmonic potentials, still axially infinite, and sectionals devoted to 
the full GPE calculation of a three dimensional elongated anharmonic trap. 



II. GROUND STATE AND BOGOLIUBOV QUASI-PARTICLES IN AN INFINITE CYLINDER 

Since the traps discussed in this paper are cigar shaped, the dominant contribution to the inhomogeneity originates 
from the trap geometry in the radial direction. Thus we first assume an infinite trap in the axial direction z. We solve 
the linearized GPE for a flat-bottom trap, uniform in the axial direction. In order to compare with our harmonic 
trap, we fix the chemical potential /i/ft. — 2.06 kHz (in the Thomas-Fermi approximation) and the number of atoms 
per unit length N = 10^ atoms per L = 52.2/xm to coincide with our BEG experimental parameters [j. 

Our flat-bottom trap potential, is a cylindrically symmetric hard-wall confinement in the radial direction: 

I oo otherwise 

where Rtf is the radial Thomas- Fermi radius, which is determined by the constraints on fi and N/L. We calculate 
the Bogoliubov spectrum for the flat-bottom trap following the derivation of for harmonic traps. Briefly, in the 
linear regime, the condensate wavefunction can be expanded as: 

*(r, t) = e''(''/'*)*{V'o(i-) + ^[c,w,(r)e-*"^* + c*v*{v)e''^^']} (2) 



where the c^ 's are constants, and u and v are the quasi-particle amplitudes, that obey the following orthogonality and 
symmetry relations: 

dr^u* = 5,; /dr{u.«,-..u,} = (3) 

The wavefunction expansion in equation |2] is inserted into the GPE: ih-^^f — — h Vt ^, where g = 

'^^^ ° is the coupling constant, m is the atomic mass and a is the s-wave scattering length. For a cylindrical condensate, 
infinite in the axial direction, we define ipo{p,z) — (j>o{p)/L, {u,v)n,k{p, z) = -^e^''^{u,v)n,k{p) as the new ground 
state wavefunction and the quasi-particle amplitudes, respectively. Gathering terms to zero order m, v yields the 
stationary GPE: 

( - ^ + Vt{p) + ^Mp)A Mp) = Mp) (4) 

\ 2m L I 

The next order gives the coupled equations: 

riUJn,kUn,k[p) = [ h h Vt(P) - —j—MPmn.kiP) + ^ 00 (/0)«n,fc (p) (5) 



, , , ^^V2 h^k^ 2gN 2, M / X .9^ 2/ ^ 

- TUVn,kVn,k(P) = [ ^ + + Vt(P) ~ P+ -J^n{p)\Vn,k{P) + 00 (P) (6) 

also known as the linearized GPE or Bogoliubov equations. 

First, we find the ground state wavefunction from equation^ by applying imaginary time evolution [loj to the 
Thomas- Fermi wavefunction ■ Figure ^ shows the ground state wavefunctions for the flat-bottom trap (solid 
line) and for an axially infinite harmonic trap (dashed line), with a confining potential of Vr = ^mw^p^ where 
Up — 27: X 226Hz 0. Indeed, the harmonic trap wavefunction is similar to the Thomas Fermi ansatz, and the flat 
trap has a ground state in which the density drops to zero at Rtf over a length scale ^ = (87rn(0)a)-i/2 

where 

n(0) is the peak density. 
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FIG. 1: Ground state radial wavefunctions of axially infinite harmonic (dashed line) and flat-bottom (solid line) traps obtained 
by solving the stationary GPE with imaginary time evolution. 





FIG. 2; Quasi-particles amplitudes u„^k and effective potentials, (a) Harmonic trap, uo,fc (solid line) and ui^k (dashed line), 
(b) Flat-bottom trap, uo,k (solid line) and ui^k (dashed line), (c) Harmonic trap, effective potential (solid line) and Bogoliubov 
frequencies (dashed line) . (d) Flat-bottom trap, effective potential (solid line) and Bogoliubov frequencies (dashed line). 



In order to find the Bogoliubov frequencies Wn,k and the quasi-particles amplitudes (u,w)n.fc, the relaxed wave- 
function is substituted in equations and El The Bogoliubov equations and El are an eigenvalue problem and are 
solved by expanding the functions u{p) and v{p) using a complete set of zero order Bessel functions. The problem is 
thus reduced to a matrix diagonalization and is numerically solved ;5,]. We take the excitation momentum k to be 
k = 8.06/ir7i~^ as in :6j. 

The calculated lowest order quasi-particle amplitudes iio,fe (zero order, solid lines) and ui^k (first order, dashed 
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lines) are shown in figure l^fa) and (b), for the harmonic and flat-bottom trap, respectively. For the n*'' order radial 
modes, both u and v have n nodes. 

Comparing figures Efa) and (b) reveals a dramatic difference. For the flat bottom trap (figure Efb)) uo.fc and ui^k 
are mainly centered around p — Q whereas for the harmonic trap (figure |2Ia)) they are significantly shifted. 

To qualitatively explain this effect we approximate equation |5l as: 



huJn,kUnAP) = [ TT-^ + 1^ ^ ^t(p) " M + '/'o (p) (7) 

Zm 2m L 

where we neglect the contribution of v, which are small for the A:^ > 1 regime considered here. The two Bogoliubov 
equations Inland El are thus reduced to a Schrodinger equation for the excitations, with the effective potential Ve// = 

^^^ + Vt(p) — /i+ ^^00 (p) ^M- figure|3c) and (d) we plot the effective potential for the harmonic and flat-bottom 
trap, respectively. The dashed lines are the lowest eigenvalues ujn,k of equations and El without neglecting v. The 
effective potential of the harmonic trap can be approximated for the lowest modes as a harmonic potential, leading 
to u{p)n,k which resemble eigenstates of the harmonic trap for u;„ sufficiently deep in the effective trap. However, 
for eigen-energies hujn,k which approach the effective trap height, there is a nonvanishing amplitude that interferes 
constructively near p — due to radial symmetry, so that the amplitude is not concentrated in the well. This happens 
even for the lowest mode in the flat trap. The well in this case is extremely narrow, originating only from the spatially 
smoothing of the ipQ at the walls, supporting no "bound states" concentrated in the well. 



III. BOGOLIUBOV BRAGG SPECTROSCOPY OF A FLAT-BOTTOM TRAP 



Once the quasi-particles amplitudes and the Bogoliubov frequencies are obtained, they can be used to derive the 
Bogoliubov Bragg spectrum, using the quasi-particle projection method |5^. Briefly, in this method the cj coefficients 
in equation 121 are taken to be time-dependent and then substituted into the GPE, which includes now also the Bragg 
external potential Vscos^kz — tot), where Vb is the strength of the Bragg potential. The momentum given to the 
condensate, in the Bragg process, defined as: 



is found to be Q: 



rj u\ NqV^t^ ^ 3 szn[(^„,fc - uj)t/2] ^ sin[{uJnM + ^^)t/2] ^ .o^ 
^z[t) — jT y vVnM [( — 7 — j - ( — ; r-p: — j 



where Wn,k = 27r / dpp[u*^ ^.{p) + ti* j,(p)]0o(p)- Wn,k can be considered as a weight function that determines the 
weight of each radial mode. This weight is given by the overlap between the ground state and the quasi-particles 
amplitudes {u,v)n,k- 

Figure |2t and show the linearized spectrum (solid lines) obtained from equation El for the k — 8.06pm~^ 
momentum, after application of a 4.2 ms Bragg pulse, in the harmonic and fiat-bottom trap, respectively. The 
harmonic trap spectrum of figure EJa) shows four resolved radial modes. It is similar to the Bragg spectrum obtained 
from full GPE calculation, for harmonically trapped elongated BEG (figure 2a in 0), as also noted in The lowest 
radial mode has a small overlap with the ground state, due to its ring shape, as seen in figure[2Ia). This ring shape can 
be explained using energy considerations. Due to exchange symmetry, the mean- field energy of an excitation is twice 
that of the ground state, making it more favorable to concentrate the quasiparticle in a ring around the ground-state. 
In the effective potential picture, the minimum of the well is near the Thomas-Fermi radius for exactly this reason. 
The spectrum may thus be qualitatively described by the following argument. For the lowest radial mode, the overlap 
of the excitation wavefunction, which is approximately gaussian, with the condensate wavefunction is small. As the 
radial number n increases, but remains small enough for the effective potential to be approximated by a harmonic 
oscillator the overlap with the ground state increases since the quasiparticle wavefunction is wider. Upon further 
increasing n, the oscillations in u and v wash out any overlap with the condensate. 

From figure EJb) it is apparent that the flat-bottom trap supports mainly the first radial mode, strongly suppress- 
ing the excitation of higher order modes, even though they are valid solutions to the Bogoliubov equations. This 
suppression is due to the small overlap Wn^k between the quasiparticle amplitude for n > and the ground state, as 
seen from in figures ^ (solid line) andj^Jb). There are no states deep in the well of the effective potential, and since 
only UQ^k is nodeless it has the dominant overlap with 00- In the Thomas- Fermi limit of an effective potential of a 
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FIG. 3: Bragg spectra for Bogoliubov (extremely weak Bragg pulses) excitations, (a) Harmonic trap, (b) Flat-bottom trap. 
Solid line corresponds with the linearized GPE solution, dashed line corresponds with the LDA. The y axis is the quantity 
Pz/{NViTik) in units of {hujp)-^. 

cylindrical well with infinite walls, the quasi-particle states are Bessel functions of order 0, Jq. This approximation 
yields a significant overlap between u and 0o for a number of radial modes. The overlap decays as l/n for large n. 
Only when considering the relaxed wavefunction, which is closer in shape to the function Jq, does the lowest radial 
mode become dominant. 

It is insightful to compare these results to the inhomogeneous line shape calculated within the LDA shown 
as dashed lines in figure 13 The LDA lineshapes were obtained by integration over the structure factor weighted 
by the inhomogeneous density distribution of the relaxed wavefunctions. The simplified LDA picture, whose basic 
assumptions are invalid here indeed fails to produce the radial mode structure. Surprisingly, it is seen from figure |31 
that it describes the envelope of the spectrum for both traps quite well, suppressing modes outside it . 

IV. OTHER CYLINDRICAL ANHARMONIC TRAPS 

So far, an idealized infinite cylindrical trap was studied. However, such a trap cannot be realized and therefore 
we consider experimentally achievable traps of finite high-order potentials . Here we show that even rather steep 
potential may yield deviations from the single-mode spectrum of the flat-bottom trap. This is done by evaluation 
of the Bogoliubov excitation spectrum for k = 8.06/im~"'^, in axially infinite traps with anharmonic radial potentials 
of the form Vt{v) — KpP^ with 2 < p < 20, using the procedure described above to solve the linearized GPE. The 
number of atoms per unit length is chosen to be the same as in section^ The chemical potential is the same for all 
traps, but different than that in section ITll 

The resulting spectra for the potential of /o^° and are shown in figures 0Ja) and (b). A strong suppression of 
modes is apparent, in comparison to the number of modes in the harmonic trap (figure |3f a)). The appearance of a 
second radial mode in the p^'^ and not in the p^^ trap seems surprising considering the flatness of both potentials. 
The qualitative difference between the two traps, which leads to such different response can be seen in the effective 
potentials plotted in figure In the effective potential picture of the p^'~' trap (figure ISJa)) we discern two trapped 
energy states, whereas for the p^*^ trap (figure E^b)) we observe only one "bound" eigen-state, while the higher 
modes are now "free". This leads to a dramatic change in the Bragg coupling matrix elements, suppressing all 
higher order radial modes, although these modes are within the allowed LDA envelope. Note that the this picture 
is only approximate due to the neglect of the v terms. Interesting corrections are expected at lower fc^, where the 
V contribution is no longer negligible. Also, the effective Schrodinger equation is radial, leading to constructive 
interference near p = 0. 
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FIG. 4: Bragg spectra for Bogoliubov excitations, (a) p^'^ trap, (b) trap. The y axis is the quantity Pz / {NVgUk) in units 

of {flUJp)~^. 




FIG. 5: Effective potentials (solid lines) and linearized GPE Bogoliubov frequencies (dashed lines) (a) trap, (b) p^° trap. 



Wc calculated the Bragg spectra for the Bogolibuov excitations for other values of and found a monotonic 
reduction in the weight of the higher order radial modes with increasing steepness of the radial potential. 



V. REALISTIC THREE-DIMENSIONAL ANHARMONIC TRAPS 



Finally, we compare between the spectra obtained above with the linearized GPE for cylindrical traps, with no axial 
confinement, and the spectra obtained for realistic elongated condensates confined in similar experimentally-realizable 
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FIG. 6: Bogoliubov spectra for the potential Vrir) = Kpp^'^ + k^z^^ , using the full GPE. 



three-dimensional anharmonic traps. We consider one specific optical trap generating a potential Vt{j^) = Kpp^^+K^z^'^ 
which is very similar to the one used experimentally for the atom-optics billiard of |l3|. 

For this we numerically solve the full GPE ihdti^ = {-?i^V2/(2m) + F+.g|V'P}'0, withy(r,t) = VTir) + VBCOs{kz- 
Lot). We evolve on a two dimensional grid (using the cylindrical symmetry of the problem) y. Nr = 2048 x 32, 
using the Crank-Nicholson differencing method. We ignore the influence of gravity to keep the cylindrical symmetry 
of the problem. Experimentally, gravitation may be compensated for using a magnetic field |l4j . 

The Bogoliubov Bragg spectrum thus calculated is shown in figure |^ Comparison between figure IHl and figure 
indicates a qualitative similarity between the spectra of the "real" three-dimensional traps and that of their 
idealized axially infinite (two-dimensional) version, even though figure includes additional inhomogeneity due to 
axial confinement, as well as possible finite-size (Doppler) effects 0,l3- Thus we see the essence of the inhomogeneous 
broadening is indeed captured by the radial dependence of the trap. 



VI. CONCLUSION 



In conclusion, we propose and study a possible approach to reduce the inhomogeneity of trapped condensates by 
changing the functional form of the trap. We show a reduction of the inhomogeneous line-shape and suppression 
of the multiple peaked structure in the Bragg Bogoliubov spectra by using anharmonic traps with steep walls. The 
suppression is found to be monotonic with the steepness of trap potential walls. Linearized calculations for simpli- 
fied two-dimensional cylindrical traps give qualitative understanding of the suppression of high-order radial modes, 
due to decreasing overlap with the ground-state wave function. Further LDA calculations we performed, suggest a 
dramatic broadening of the spectra in the presence of gravity, emphasizing the necessity to compensate for gravity in 
experimental realization of the anharmonic traps. Finally, full GPE simulation of the spectra with previously used 
optical traps confirm the experimental feasibility of this approach. The reduction in inhomogeneous broadening opens 
the possibility to study the more intrinsic homogenous broadening mechanisms. Experimental realization of BEC 
confined by anharmonic traps should also yield new experimental opportunities, such as the sharpening of the phase 
transition from superfluid to Mott-insulator and rich variety of vortex phases [l6j . 

This work was supported in part by the Israel Ministry of Science and the Israel Science Foundation. The authors 
wish to express their gratitude to C. Tozzo and F. Dalfovo for their assistance. 
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